Off-Fermi Shell Nucleons in Superdense Asymmetric Nuclear Matter 
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Recent observations of the strong dominance of proton-neutron (pn) relative to pp and nn short- 
range correlations (SRCs) in nuclei indicate on possibility of unique new condition for asymmetric 
high density nuclear matter, in which the pp and nn interactions are suppressed while the pn in- 
teractions are enhanced due to tensor interaction. We demonstrate that for sufficiently asymmetric 
case and high densities the momentum distribution of the smaller p-component is strongly deformed 
with protons increasingly populating the momentum states beyond the Fermi surface. This re- 
sult is obtained by extracting the probabilities of two-nucleon (2N) SRCs from the analysis of the 
experimental data on high momentum transfer inclusive electro-nuclear reactions. We fitted the 
extracted probabilities as a function of nuclear density and asymmetry and used the fit to estimate 
the fractions of the off-Fermi shell nucleons in the superdense nuclear matter relevant to neutron 
stars. Our results indicate that starting at three nuclear saturation densities the protons with frac- 
tional densities i will populate mostly the high momentum tail of the momentum distribution while 
only 2% of the neutrons will do so. We discuss the implications of this condition for neutron stars 
and emphasize that it may be characteristic to any asymmetric two component Fermi system with 
suppressed central and enhanced short-range tensor interactions between the two components. 
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The recent experiments on high-momentum transfer 
semiexclusive reactions 0,0 in which the struck nucleon 
from the nucleus is detected in coincidence with the re- 
coil nucleon from SRC, found striking disbalance be- 
tween pn and pp/nn correlations. They found that pro- 
tons struck from the nucleus with initial momenta of 
kp < p < 600 MeV/c, in the 92% of the time emerge 
from the pn SRC, while the pp and nn SRCs are sig- 
nificantly suppressed, contributing only « 4% to the 
high momentum part of the nucleon momentum distri- 
bution in nuclei. This disbalance was understood based 
on the dominance of the NN tensor interaction in the 
300 < p < 600 MeVVc momentum range relevant to 2N 
SRCs in nuclei [l|, [3, The tensor- interaction dominance 
is due to the fact that in this momentum range corre- 
sponding to inter-nucleon distances of 1 Fm the NN 
central potential is crossing the zero due to transition 
from attractive to the repulse interactions. As a result 
at these distances overall NN potential is dominated by 
tensor interaction, which results to the suppression of s- 
channel isostriplet pp and nn interactions and enhance- 
ment of the interaction in isosinglet, L—2, pn channel. 
The resulting picture for the nuclear matter consisting of 
protons and neutrons at densities in which inter-nucleon 
distances are ~ 1 Fm is rather unique: it represents a 
system with suppressed pp and nn but enhanced pn in- 
teractions. 

The goal of the present study is to understand how 
the high momentum part of the momentum distributions 
of protons and neutrons are defined in the high density 
nuclear matter under the above described conditions. 

New Relation between High Momentum p- and 
n-distributions in Nuclei: Due to short range na- 
ture of NN interaction the nuclear momentum distribu- 
tion, n^{p), for momenta, p, exceeding the characteristic 
nuclear Fermi momentum kp is predominantly defined 



by the momentum distribution in the SRCs. There is a 
rather large experimental body of information indicating 
that for the range of fc^ < P < 600 MeV/c the SRCs are 
dominated by 2N correlations, which consist of mainly 
the pn pairs (for recent reviews see [ll|, [l^] ) ■ 

In recent work[l^ based on the dominance of the pn 
SRCs we predicted two new properties for the nuclear 
momentum distributions a.t kp < p < 600: (i) There is 
an approximate equality of p- and n- momentum distri- 
butions weighted by their relative fractions in the nucleus 
Xp and Xji . 

Xp-ripip) ^ x„n^{p), (1) 
with J n^^^{p)(Pp = 1. (ii) The probability of proton 
or neutron being in high momentum NN SRC is inverse 
proportional to their relative fractions and can be related 
to the momentum distribution in the deuteron nd(j)) as: 
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-a2{A,y) ■ nd{p), 
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where 02 (A, y) is interpreted as a per nucleon probability 
of finding 2N SRC in the given A nucleus 0-01 and the 
nuclear asymmetry parameter is defined as y = 1 1 — 2xp \ . 

The above two properties are obtained assuming no 
contributions from pp, nn as well as higher order SRCs. 
They follow from the assumption that the whole strength 
of nuclear high momentum distribution as well as per 
nucleon probability of proton and neutron to be in the 
SRC is defined by the same pn correlation. In Ref. (l3| we 
demonstrated that these properties are seen in the direct 
calculations using realistic ^He wave function. 

Since the SRC is defined by local properties of nuclei, 
one expects that the A dependence of 02, in Eq.([2|) is re- 
lated to the nuclear density, i.e. a2{A, y) = 02 (p, y). This 
could allow us to estimate the high momentum part of 
the nucleon momentum distribution not only for finite [l3l| 
but also for infinite nuclear matter. 
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The Parameter a2 and A{ee')X Processes: In prin- 
ciple, any nuclear process which probes high momentum 
nucleons in nuclei should allow an extraction of 02 {A, y) . 
One of such processes are high momentum transfer inclu- 
sive A{e, e')X reactions measured in special kinematics 
in which electron scatters off a deeply bound nucleon 
having large momentum in the nucleus. Two param- 
eters, 4-momentum transfer square —Q^ and Bjorken 
XB-i = — allow us to select these kinematics. In- 
troducing the parameter a, which defines (A times) the 
light cone momentum fraction of nucleus carried by the 
interacting nucleon, within impulse approximation: 



XBj 
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where 4-momentum of the virtual photon is defined as 
(90, 93) with q± = go ± 93- Also pi+ and rh represent the 
-component and mass of the bound initial nucleon. 
From Eq.© in the limit of > m% such that ^ < 1 
the condition: a « xsj is satisfied and choosing xbj > ^ 
will select a nucleon in the nucleus that carries momen- 
tum fraction more than that of the stationary nucleon. 
This observation is the basis of the 2N SRC model [H-Q, 
according to which at ^ nij^ and 1.4 — 1.5 < xbj < 2 
the interacting nucleon needs to acquire a substantial mo- 
mentum fraction from the nucleon with which it is in a 
short-range space-time correlation. The expectation that 
the 7* A*" interaction in SRC will be weakly influenced by 
the long-range mean- field of A-2 residual nucleus resulted 
to the prediction of the onset of plateau in the ratios of 
A(e,e')X and d(eje')X cross sections at xbj > 1-4 — 1.5 
and » w^j^-Q- First, such plateau was observed in 
Ref. T| and later was confirmed in new experiments [8l-[l0| 
for wide range of nuclei. The measurements also con- 
firmed that the onset of the plateau depends on Q"^ and 
sets in at > 1.5 GeV^ as it was predicted in the 2N 
SRC model (see e.g. [1, i, [H Hj]). It is worth noting 
that models in which the x > 1 cross section is attributed 
mainly to the final state interaction of the struck nucleon 
with the residual nucleons is in disagreement with the ob- 
served plateau and its onset being a function of (for 
detailed discussion of FSI effects see Ref.[l2[). 

Based on the expectation that A{e, e')X probers 2N 
SRCs, from Eq.([5]) one observes that 

2(7,, A ■ , dcr 
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with 
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for those values of xbj and that the measured ratio 
of the cross sections exhibits the plateau. 
Extraction of 02 (A): We analyzed the compilation 
of the world data on inclusive A{e, e')X reactions from 
Ref.P,!!!]. Only the data for d, ^He, ^He, ^Be, ^^C, 
^''AZ, ^^Fe, ^'^Cu and ^^^Am nuclei satisfied the criteria of 
XBj > 1-5 and > 1.5GeV^. We first constructed the 
data matrix for central and x spanning the following 
values: = 1.75, 2.25, 2.75 3.25 and xbj = 1.55, 1.65, 
1.75. For each pairs of Q^, XBj we averaged the CeA cross 
sections and their errors with AQ^ — ±0.25 GeV^ and 



TABLE I: The results for a2{A,y) 



A 


y 


This Work Ref.f7] 


Ref.f8, 9] 


Ref.|40l 


^He 


0.33 


2.07±0.08 


1.7±0.3 




2.13±0.04 


^He 





3.51±0.03 


3.3±0.5 


3.38±0.2 


3.60±0.10 


^Be 


0.11 


3.92±0.03 






3.91±0.12 







4.19±0.02 


5.0±0.5 


4.32±0.4 


4.75±0.16 


^^Al 


0.037 


4.50±0.12 


5.3±0.6 






56 pg 


0.071 


4.95±0.07 


5.6±0.9 


4.99±0.5 






0.094 


5.02±0.04 






5.21±0.20 


"^Au 


0.198 


4.56±0.03 


4.8±0.7 




5.16±0.22 



/S.XBj — ±0.05. Along with the averaged cross sections 
we estimated the average values of relevant kinematic 
variables, 7 for each bin according to: 



The 02 is estimated for each 
a2{A,y 



XBj 



<7eA{xA,Q 



<7edixd, Qi, Oe.d) 



bin as: 



i?, 



(5) 



(6) 



where a;A(d), Q'A(d) ^^'^ Oe,A(d) are average values for nu- 
clei and d defined according to Eq. ([S]) . The factor R uses 
the theoretical calculation of d{e, e')X reaction [l^[l3| to 
correct for the misalignment of averaged x, and for 
nuclei A and d in the following form 



R 



aZ{Qlxd,e,,d) 
af^{Q\,XA,e,,A) 



(7) 



where cr*'* is a model calculation of the cross sections. 

The results of 02 for the above mentioned nuclei 
are given in Fig. 1 and in Table I together with the 
previous0-[ll and recent[l3| estimates. Our results are 
somewhat lower than that of Refs.[7-'To'], and together 
with Ref.[l3] they agree with the earlier indication^] that 
02 decreases for heaviest nuclei due to larger asymmetry 
y which is in agreement with the observation of the sup- 
pression of nn/pp vs pn SRCs. 

Fitting of 02(^.5 y): We now use the extracted values 
of 02 to fit them in the parametric form: 



a2(A,2/) = a2(A0)/(y). 



(8) 



The justification for the factorization of A and y depen- 
dences follows from the fact that the asymmetry depen- 
dence of 02 is due to its proportionality to the number of 
the pn pairs per nucleon. Thus one expects same function 
f{y) for nuclei with different A. 

First, we fit a2(A, 0)s for symmetric nuclei. However 
we have only two data points for 02 (A, 0): ^TJe and "^^C. 
To be able to fit the <i2(A, 0]s for the range of A > 12 we 
use the approximation [Sl la. Ill | : 



a2(A,0) = C / p\{r)d\ 



(9) 
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FIG. 1: (a): The A dependence of 02. The triangles are the 
extracted values of 02, squares - scaled {p\) (see the text), 
dashed line - 02(^,0), solid line - final fit for 02 (A, j/). (b) 
The y dependence of the asymmetry function. 



where PAif) is the nuclear matter density with 
J PAir)cfir — 1. This relation follows frora the propor- 
tionality of 2N SRCs of finding two nucleons at the same 
position which is related to the second order of the nu- 
clear matter density function. The fact that the matter 
density depends weakly on nuclear asymmetry at large 
A justifies the use of this ansatz for fitting 02 (A, 0). 

We calculated the {p\) = J p\{r)(Pr for differ- 
ent nuclei using the pA parameterizations extracted 
from the experimental measurements of nuclear charge 
densities [l^. Normalizing the calculation of (p^) for ^^C 
to the extracted 02(12) from Table I yields. 

C = 49.1 ±2.6. (10) 
The 02 (A, 0) "data", generated in this way are presented 
in Fig.l. Combining these a2(A,0)s with the extracted 
flaC^iJe) (Table I), we then fit a2(A,0)'s for the whole 
range of A as it is presented in Figllja) (dashed line). 

Next, we use the 02(^,0) fit and Eq.([5]) to extract the 
six data points for f{y) from the measured 02 's for asym- 
metric nuclei (FiglTJb)). To be able to make a meaningful 
fit we supplement these points with the following condi- 
tions that have transparent physics interpretations: (a) 
positiveness condition, /(y) > for any y; (b) bound- 
ary condition, /(O) = 1 and /(I) = 0, where the sec- 
ond condition follows from the approximation in which 
we neglected pp and nn SRCs; (c) since the Ay 
limit is valid only in the case of A — > 00 one obtains 
that /'(O) = /'(I) = 0. The general ansaz for f{y) with 
minimal number of free parameter satisfying conditions 
(a) — (c) can be presented in the following form 

f{y) = (1 + (6 - 3)y2 + 2(1 - h)y^ + by^)F{y), (11) 
where the additional correction function F{y) accounts 
for the non-smoothness of the asymmetry curve at y < 
0.15. The best fit is obtained for 6 « 3. Combining 
02(^,0) and f{y) fits in Eq.(|Sl) we obtain the final fit 
which is the solid line in FiglIJa). 



Extrapolation to Infinite and Superdense Nuclear 
Matter: The obtained fit in Eq.([8l) allows us to estimate 
02(^1 y) for infinite nuclear matter since Eq.(|ni) converges 
at A — > cx) and /(y) is finite by definition. The estimate 
for the symmetric nuclear matter at saturation densities 
Po can be obtained using the relation between the nuclear 
radius and A, R — r^) ■ A^ , which yields 



{P 



2\INM 
I sym 



PA,symir)d^r 



47r 



plr^ « 1.4 fm \ 



(12) 

where we use po = 1-6 fm and = 1.1 fm. From 
Eqs.([S]) and (flUl) we obtain for symmetric nuclear matter 
at saturation density: 



a2(po,0) « 7.03 ±0.41, 



(13) 



which is quantitatively in agreement with the 02 esti- 
mated from the y scaling analysis of the A(e,e')X data ex- 
trapolated to infinite nuclear matter [23| which yields [2l| 
02 « 8.0 ± 1.24. Note that our estimate gives the lower 
limit for 02 due to the neglection of the nn and pp SRCs. 




neutrons 



P/Po 



FIG. 2: Density dependence of the fraction of off-Fermi-shell 
nucleons in a^r, = ^ matter. 



Next we consider asymmetric nuclear matter. We com- 
bine Eqs. (|T2|) and ([TT|) into Eq.(l8l) to estimate 02 (p, y) for 
given values of nuclear density and asymmetry y. As an 
example of the application of 02 {p, y) , we estimate the 
fraction Pp/n of off-Fermi-shell nucleons in the /3 equi- 
librium e — p — n superdense asymmetric nuclear matter 
using the relation (see Eq.(l2|): 



Pp/n{A, y) 



1 



2x. 



p/n 



a2{A,y) nd{p)d'p. (14) 



For the asymmetry y, in our estimates we use the thresh- 
old value of Xp = i (corresponding to y = |) below of 
which the direct URCA processes: 

n — !■ p + + p + e^^n + Ve (15) 

will stop in the standard model of superdense nuclear 
matter consisting of degenerate protons and neutrons [22I]. 
Estimating the Fermi momenta of protons and neutrons 
in Eq. ([T4l) with kp.N = {iT^"^ x n p)^ , in Fig. 2 we present 
the off-Fcrmi-shell fractions of protons and neutrons as a 
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function of nuclear density. The most interesting result 
of these estimates is that in equilibrium pn SRCs move 
the large fraction of protons above the Fermi-shell: at 
Spo densities half of the protons will be ofF-Fermi-shell 
while at p > 4.5/9o all the protons will populate the high 
momentum tail of the momentum distribution. The situ- 
ation however is not as dramatic for neutrons, with only 
few% of neutrons populating the high momentum part 
of the momentum distribution. 

Possible Implications for Nuclei and Neutron 
Stars: Our main observation is that with an increase of 
nuclear asymmetry the lesser component become more 
energetic. This is confirmed 13.] by direct calculation 
of the average kinetic energies of proton and neutron 
using realistic wave funcion of ^He, in which case one 
expects neutrons to be more energetic than protons: 
(r„) = 18.4 MeV and (Tp) = 13.7 MeV. For nuclei with 
large A(> 40) one expects protons to be more energetic 
than neutrons, with larger fraction of protons occupy- 
ing high momentum tail of the momentum distributions. 
This may have several verifiable implications for large A 
nuclear phenomena [isj . 

Our observation may have more dramatic implications 
for the dynamics of neutron stars. Some of them are: 

- Cooling of a Neutron Star: Large concentration of pro- 
tons above the Fermi momentum will allow the condition 
for Direct URC A processes Pp+ Pe > Pn [22] to be satis- 
fied even ii Xp < i. This will allow a situation in which 
intensive cooling of the neutron stars continues well be- 
yond the critical point Xp = ^ (see also Rcf. 11]). 

- Superfluidity of Protons: Transition of protons to the 
high momentum tail will smear out the energy gap which 
will remove the superfluidity condition for the protons. 

- Protons in the Neutron Star Cores: The concentration 
of protons in the high momentum tail will result in proton 
densities Pp ^ Pp ^ k^^p- This will favor an equilibrium 
condition with "neutron skin" effect in which large con- 
centration of protons populates the core rather than the 
crust of the neutron star. This and the proton superflu- 
idity condition violation may provide different dynamical 
picture for generation of magnetic fields in the stars. 

- Isospin locking and the stiff equation of state of the 
neutron stars: With an increase in density more and 
more protons move to the high momentum tail where 
they are in short range tensor correlations with neutrons. 
In this case one would expect that high density nuclear 
matter to be dominated by configurations with quantum 



numbers of tensor correlations (S* = 1,7 = 0). In such 
a scenario protons and neutrons at large densities will 
be locked in the NN iso-singlet state. This will double 
the threshold of inelastic excitation from NN NA to 
NN — > AA{NN*) transition thereby stiffening the equa- 
tion of state which is favored by the recent large neutron 
star mass observation [235. 

Possible Universality of the Obtained Result: 
Our observation is relevant to any asymmetric two- 
component Fermi system in which the interaction within 
each component is suppressed while the mutual interac- 
tion between two components is enhanced. It is inter- 
esting that the similar situation is realized for two-fermi- 
component ultra-cold atomic systems 24] but with the 
mutual s-statc interaction [2 7j. One of the most intrigu- 
ing aspects of such systems is that in the asymmetric 
limit they exhibit very rich phase structure with indica- 
tion of the strong modification of the small component of 
the mixture psl [26| . In this respect our case is similar to 
that of ultra-cold atomic systems with the difference that 
the interaction between components has a tensor nature. 

Limitations and Outlook: Our analysis has several 
limitations: One is that we neglected the contributions 
from isotriplet 2N as well as 3N SRCs. Even though 
the statistical errors in the extraction of 02 are small 
(Table I) , an additional errors are accumulated due to the 
fitting procedure, especially for the asymmetry function 
f{y)- We estimate the overall error in the extrapolation 
procedure at ~ 30%. 

Finally, the procedure of extraction and fitting of 
0-2 {p, y) can be significantly improved with the new high 
and XBj > 1 experiments covering widest possible 
range of A and y. The semiinclusive A{e, e'NN)X data 
will allow an inclusion into the analysis a contribution 
from pp and nn SRCs. Measurements at xbj > 2 do- 
main will allow also to obtain similar estimates for 3N 
SRCs. Inclusion of all these effects into the analysis will 
further increase the magnitude of high momentum frac- 
tion of the protons. Thus our present results represent 
most probably the lower limit of the fraction of off- Fermi 
shell protons in high density nuclear matter. 
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